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ABSTRACT 


Context. Asteroseismology provides a powerful way to constrain stellar parameters. Solar-like oscillations have been observed on 
subgiant stars with the Kepler mission. The continuous and high-precision time series enables us to carry out a detailed asteroseismic 
study for these stars. 

Aims. We carry out data processing of two subgiants of spectral type G: KIC 6442183 and KIC 11137075 observed with the Kepler 
mission, and perform seismic analysis for the two evolved stars. 

Methods. We estimate the values of global asteroseismic parameters: Av = 64.9 ± 0.2 pHz and v^ax = 1225 ± 17 pHz for KIC 
6442183, Av = 65.5 ± 0.2 pHz and Vn,ax = 1171 ± 8 pHz for KIC 11137075, respectively. In addition, we extract the individual mode 
frequencies of the two stars. We compare stellar models and observations, including mode frequencies and mode inertias. The mode 
inertias of mixed modes, which are sensitive to the stellar interior, are used to constrain stellar models. We define a quantity dv^_p 
that measures the difference between the mixed modes and the expected pure pressure modes, which is related to the inertia ratio of 
mixed modes to radial modes. 

Results. Asteroseismic together with spectroscopic constraints provide the estimations of the stellar parameters: M = l.OTj^gQ^Afo, 
R = 1.66+p/?o and t = 8.65!^;*^ Gyr for KIC 6442183, and M = 1.00!“;°[Mo, R = 1.63+°;°[7?o and t = 10.36+°;°^ Gyr for KIC 
11137075. Either mode inertias or could be used to constrain stellar models. 
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1. Introduction 


Asteroseismology provides a useful tool to probe stellar inte¬ 
riors, test internal phy sical processes and obtain stellar 


eters accurate! 


iKallinger et al 


•nal phy sical processes and obtain stellar param- 
[ y (e.g . lEggenbergeret alJ 2004I; iBi et al.l [20081: 
T 120101 : iMontalban et alJ 120131) . Many solar-like 


stars have been observed continuousl y and precisely with space 
missions such as C oRoT and Kepler (lAppourchaux et al.ll2008t 
iBorucki et alJl200^ . ushering in a golden age for the asteroseis¬ 
mic analysis of oscillating stars. 


The oscillations of evolved stars include mixed modes, 
which behave as pressure modes (p-modes) i n the enve¬ 
lope and gravity mod es (g-modes) in the core (lOsa El [ml 
lAizenman et al.l[l977t) . Mixed modes have been used to distin¬ 
guish between red clump stars and RGB stars dBedding et alJ 
1201 It iMosser et al.ll201 fa) , and to monitor stellar evolution sta- 
tus from the main sequence to the asymptotic giant branch 
dMosser et alJl201^ . The signature of avoided crossings raises 
the exciting possibility that detailed modeling of the star will 


provide a very precise determination of its age 

Gilliland et a! 

2010al:IChaolin et al.ll2010l:lMetcalfe et al.ll2010l 

Benomar et a! 

2012HBeddingll20l4. 


As discussed 
Christensen-DalsgaardI 


by iDziembowski et al 


d2003) 


( 2001 ); 

d2009h : 


_ _ iDupret et a!. _ 

Benomar et al.l d2014l) : iDatta et ^ d2015[) . the g-dominated 


non-radial modes have much higher inertia than the radial 


modes, whereas the p-dominated modes do not. Because the 
mixed modes possess the properties of partial g-modes, the 
inertia of mixed modes is usually higher than that of p-modes. 
The inertia of mixed modes could provide powerful constraints 
on stellar models. 


Subgiants are the bridge between the main sequence (MS) 
and red giants on the H-R diagram. The study of subgiants 
could provide valuable insight into stellar structure and evolu¬ 
tion. When stars leave the main sequence stage, the g-mode and 
p-mode frequencies overlap, which results in the appearances of 
mixed modes in post-main-sequence stars. KIC 6442183 (aka. 
‘Dougal’) and KIC 11137075 (aka. ‘Zebedee’) are at the begin¬ 
ning of the subgiant stage. The magnitudes of the stars are 8.52 
and 10.86 (KIC magnitude), which are bright enough to allow 
seismic observations. They have been observed by the Kepler 
mission continuously with a short cadence of 58.84 seconds 
(iGilliland et a! 2010bh fo r quart ers 6.1-17.2 and 7.1-11.3, re¬ 
spectively. Benomar et al.l (120131 l2014t) extracted the mode fre¬ 
quencies and measured the mode inertia ratio for KIC 6442183. 
They estimated the mass M — 0.94Mo and the radius R - 
1.60Rn for KIC 6442183 using scaling relations (iBrown et al.l 
Il991l : [^eldsen & Bedding|ll995h . With oscillation observations 
from K eple r and spectroscopic obse r vation s from iBruntt et al.l 
(1201 2h and iMolenda-Zakowicz et al.l (l2013l) . we can now con¬ 
duct asteroseismic analyses and constrain stellar parameters for 
the two stars accurately. 
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Table 1. Basic parameters of KIC 6442183 and KIC 11137075 from observations. 



KIC 6442183 

KIC 11137075 

Teff (K) 

5740 ± 70‘ 

' 5738 ±62" 

5590 

±70" 5610 ±71" 

logg (dex) 

4.03 ± 0.03' 

' 4.14 ±0.10* 

4.01 ± 

0.03" 4.10 ±0.12* 

[Fe/H] 

-0.11 ±0.06' 

' -0.12 ±0.05* 

-0.06 ± 

0.06" -0.06 ±0.06* 

KIC magnitude 

8 

.52' 


10.86' 

Av (//Hz) 

65.07 ±0.09' 64.9 ±0.2'' 


65.5 ± 0.2'' 

Vmax (//Hz) 

1160 ±4' 

1225 ± 17'' 


1171 ±8'' 

M(Mo) 

0.94' 

1 r)4+0.01rf 


^‘^'^-0.01 

R(Ro) 

1.60' 



1 ^ 


Notes. [Bruntt et^ ( 1201 2h iMolenda-Zakowicz et al.l ( 1201 3h iBenomar et all ( l2013h . who stated that typical uncertainties on mass and 
radius were a few percent This work Kepler Asteroseismic Science Operations Center: http://kasoc.phys.au.dk/ 


In Section 2, we review the recent observations and extract 
the mode frequencies for the two stars. In Section 3, we con¬ 
struct stellar models and correct the near-surface term for the 
theoretical frequencies of the models. Seismic analyses and stel¬ 
lar parameters determinations are shown in Section 4. Finally, 
the discussions and conclusions are presented in Section 5. 


2. Observations and Data Processing 


Spectroscopic observations fr om iBruntt et~^ (1201 2ll and 
iMolenda-Zakowicz et al.l (1201 3l) provided the values of metal- 
licity ([Fe/H]), effective temperature {T^s) and gravity (logg) 
for these two stars. Results from the two groups show good con¬ 
sistency, and yield preliminary constraints on the stars. In this 
paper we adopt the parameters obtained by iBruntt et'aP (12012l) 
because they took the asteroseismic logg into account during the 
spectroscopic analysis. 

We obtained the power spect r a by applying the Lomb- 
Scargle Periodogram (lLombl[T976t IScairgld 11982 ) to the short 
cadence time series corrected by KASC Worki ng Group 1 

1 WG# 1; ‘solar-like oscillating stars’) following iGarcfa et al.l 
201 ih . which is ava ilable to the Kepler A steroseismic Science 
Consortium (KASC jKieldsen et al.ll201^ through the KASOC 
database B 

The autocorrelation of the power spectrum can provide 
the estimate of perio dic informat i on, s uch as mean large 
frequency separation (iBarban et al.l l2009h . In this work, we 
estimated the large frequency separations Av and the fre¬ 
quency of m aximum power Vmax with the autocorrelation 
function (e.g. Ro xburgh & Voronts^ 120061: iRoxburghI l2009t 
iMosser & Anpourchauxl 2009h and collapsed a utocorrelation 
function (e.g. lHuber et al.ll200^ iTian et al. l2014h . These global 
oscillation parameters were evaluated as: Av = 64.9 + 0.2 yuHz 
and Vmax = 1225 + 17 //Hz for KIC 6442183, Av = 65.5 + 0.2 
//Hz and v^ax = 1171 + 8 //Hz for KIC 11137075. The global os¬ 
cillation parameters (Av and v^ax) and atmospheric constraints 
for the two stars are listed in Tabled!] 

As discussed by iMathur et al.l (1201 Ih . several methods have 
been developed to extract individual frequencies, and the key 
point is to fit a sum of Lorentzian profiles describin g all oscil¬ 
lation mod es. We followed the method de scribed by Basu et aP 
(I 2 OOOI) and iRegulo & Roca CortesI (l2007h . where each oscilla¬ 
tion mode in the power spectrum was fitted to a Lor entzian pro¬ 
file using r obust non-linear least squares metho d ( Markwa^ 
1200911200 . As discussed bv lAnpourchaux et al.l(ll998allbh . the 


* Kepler Asteroseismic 
http://kasoc.phys.au.dk/ 


Science Operations 


Center: 


noise statistic of each bin of the power spectrum has a dis¬ 
tribution with 2 degrees of freedom. However the non-linear 
least squares method assumes that this statistics is Gaussian. The 
mode extraction method adopted in this work might introduce a 
certain bias of amplitudes and linewidths, and fortunately does 
not affect the frequencies. To estimate the errors, we measured 
the individual frequencies quarter by quarter, and estimated the 
frequency and error of each mode by combining each subset. 
Asymptotic formulae show that the oscillati on mode linewidthds 
are related to Tj-ff of stars with a power law (IChanlin et al.ll200^ 
lAppourchaux et af] l2012h . The smallest linewidth of the de- 
tected modes for KIC 6442183 is of approximately 0.25 //Hz 
(IBenomar et al .11201^ . Because the difference of between 
the two subgiants is about 150 K, they will obtain the linewidths 
of the same magnitude. Therefore, the frequency resolution of a 
Kepler quarter (0.1286 //Hz) is sufficient to resolve the oscilla¬ 
tion modes for the two stars. 

This procedure provided us with independently obtained fre¬ 
quency sets with symmetric errors. We obtained 37 oscillation 
modes of degree 1-0-3 for KIC 6442183 and 26 oscillation 
modes of degree 1-0-2 for KIC 11137075. The individual 
frequencies for the two stars are listed in Tables |2] and [3] The 
oscillati on frequencies for KIC 6442183 in this work and those 
given bv iBenomar et alJ (12013h are consistent at the 1 - cr level. 
The mode frequencies of the power spectra are over-plotted on 
the echelle diagram in Figs. [3] and 0] The radial modes follow 
a vertical ridge in the figures, and the non-radial modes show 
avoided crossings, especially for the / = 1 modes. In addition to 
the 1-0-2 modes, there are five I = 3 modes identified for 
KIC 6442183, which will provide extra constraints on theoreti¬ 
cal models. 


3. Stellar Models 


To estimate the parameters of the two stars, a grid of stellar evo¬ 
lutionary mod els was constructed wit h the Yale stellar evolution 
code (YREC7: lDemaraue et al.ll2008l). We used the OPAL opac¬ 
ity table GN 93 (llglesias & Rogers 1996h . the low-temperature 
table AGS05 (iFerguson et al.l2 005l). the OPA L equation-of-state 
tables EOS2005 ( Ro gers & jjavfonovl2002l) and the Bahcall nu¬ 
clear rates (iBahcall et al. 1995 ) for microphysics. We chose the 
Eddington grey atmosphere T - t relation. We a dopted the stan¬ 
dard mixing-length theory (iBohm- Vitensei 1 9^ and overshoot¬ 
ing to treat convection. T he coefficien t of he lium and heavy el¬ 
ements diffusion is from iThoul et ^ (Il994l) . We did not take 
rotation or magnetic field into consideration in our calculations. 
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Fig. 1. Power spectrum for KIC 6442183. The black and red lines denote the power spectrum before and after smoothing to 2 yuHz. 
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Fig. 2. Same as Fig.[I]but for KIC 11137075. 


The ratio of heavy-elements to hydrogen as a mass fraction 
was estimated through the formula: 


[Fe/Hl=l„g(|)-log(|)^, 


( 1 ) 


where we adopt the value of = 0.0245 (llglesias & RogersI 
Il996h and the values of [Fe/H] from iBruntt et aP (1201 2h for 
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Table 2. Observed frequencies for KIC 6442183. Table 3. Observed frequencies for KIC 11137075. 


/ 

Frequencies (pHz) 

error (pHz) 

0 

809.63 

0.11 

0 

874.08 

0.25 

0 

937.21 

0.10 

0 

1001.04 

0.12 

0 

1065.65 

0.14 

0 

1130.50 

0.09 

0 

1195.40 

0.03 

0 

1260.32 

0.12 

0 

1325.73 

0.19 

0 

1391.69 

0.55 

1 

783.35 

0.31 

1 

840.62 

0.11 

1 

901.27 

0.14 

1 

960.26 

0.09 

1 

1002.91 

0.21 

1 

1037.60 

0.09 

1 

1096.83 

0.25 

1 

1159.96 

0.06 

1 

1224.09 

0.13 

1 

1288.75 

0.22 

1 

1354.10 

0.26 

1 

1419.52 

0.37 

2 

801.63 

0.35 

2 

868.11 

0.17 

2 

931.16 

0.30 

2 

996.92 

0.15 

2 

1059.95 

0.17 

2 

1124.73 

0.23 

2 

1189.88 

0.15 

2 

1254.61 

0.14 

2 

1321.26 

0.31 

2 

1386.81 

0.41 

3 

1150.13 

0.25 

3 

1215.93 

0.31 

3 

1280.49 

0.31 

3 

1347.26 

0.37 

3 

1412.36 

0.59 


the two stars. The values of the ratio of the heavy element 
abundance to hydrogen abundance are in the range of 
0.0165-0.0218 and 0.0186-0.0245 for KIC 6442183 and KIC 
11137075, respectively. In the model calculation, we chose the 
initia l helium abundance as Y , - 0.245+ 1.54Z, fe.g.. lDotter et alJ 
120081: iThompson et al.ll2014i) . in terms of the initial metal abun¬ 
dance. For a given mass, stellar evolutionary models depend on 
three free parameters: initial chemical compositions, the mixing- 
length parameter a and the overshooting parameter ao,,. The in¬ 
put parameters for model calculations are listed in Table-Jl] 

We calculated three groups of tracks for these two stars: (i) 
models without diffusion and without overshooting, (ii) models 
with diffusion and without overshooting and (iii) models with 
diffusion and overshooting (ao,, = 0.2). Because the observed lu¬ 
minosities of the two stars were not available, we used the large 
frequency separation Av to constrain stellar models in the H-R 
diagram. Considering the difference between Av of the models 
from the scaling relation and pulsation code, we assigned the 
large separation an error of approximately 2 fjHz in the H-R dia¬ 
gram. There are 256 and 265 tracks falling into their error boxes 
for KIC 6442183 and KIC 11137075, respectively, as shown in 
Fig.|5] 


1 

Frequencies (pHz) 

error (pHz) 

0 

949.02 

0.22 

0 

1013.37 

0.22 

0 

1078.56 

0.09 

0 

1144.34 

0.10 

0 

1209.66 

0.13 

0 

1275.11 

0.08 

0 

1341.11 

0.25 

0 

1407.67 

0.45 

1 

916.94 

0.33 

1 

978.01 

0.11 

1 

1040.74 

0.10 

1 

1101.50 

0.10 

1 

1142.65 

0.12 

1 

1180.26 

0.09 

1 

1240.80 

0.07 

1 

1304.94 

0.19 

1 

1370.07 

0.62 

1 

1436.91 

1.05 

2 

943.76 

0.47 

2 

1007.47 

0.21 

2 

1072.59 

0.21 

2 

1132.52 

0.17 

2 

1141.38 

0.23 

2 

1204.19 

0.07 

2 

1269.64 

0.20 

2 

1335.61 

0.73 



Fig. 3. Echelle Diagram with the identified modes for KIC 
6442183. The squares denote the / = 2 mode, while the circles 
/ = 0, the triangle I - 1 and the five-pointed star I - 3 modes. 


4. Asteroseismic Diagnostics 

4.1. Model calibration 

For stellar models along the tracks falling inside the error box, 
we computed smaller and more finely sampled grids around 
these models, and calculat ed theoretical m ode frequencies with 
Guenther’s pulsation code (lGuentheJl994l) . It is well known that 
there is a systematic offset between observed and computed fre¬ 
quencies that arises from improper modeli ng of the near-surface 
layers for both the Sun and solar-type stars dKieldsen et al.l20()M 
and references therein). To corr ect the near-sur f ace ter m, we fol¬ 
lowed the method described bv iBrandao et alJ d201 ll) . who cor- 
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Frequency modulo 65.5 /xHz 


Fig. 4. Echelle Diagram with the identified modes for KIC 
11137075. The squares denote the I = 2 modes, while the cir¬ 
cles / = 0, the triangle I - 1 modes. 


Table 4. Input parameters for model calculations. 


Variable 

Range 

6 

Mass (Mo) 

0.90-1.08 

0.01 

z, 

0.012-0.018 

0.001 

a 

1.70-1.90 

0.2 

^ov 

0.0-0.2 

0.2 


3-78 3.77 3-76 3.75 3.77 3.76 3.75 3.74 3.73 

logT,,, (K) logT,,, (K) 



Fig. 5. Tracks for KIC 6442183 falling into the error box in the 
H-R diagram. The left and right panels denote models for KIC 
6442183 and KIC 11137075, respectively. Red and blue colors 
in each panel denote models with mixing-length parameter a - 
1.70 and a — 1.90, respectively. 


rected the n ear-surface term for r adial modes with the method 
proposed bv iKieldsen et al.l (l2008h by fitting a power law; 


7 AN 7 AN /Fobs(n,0)'\ 

Vobs(n, 0) - Vbest(n, 0) = fl - . 

\ Fo ) 


( 2 ) 


For the mixed modes, the corrected frequencies Vcorr(n, 0 could 
be calculated through; 


Vcorr(n, /) = Vbe.st(n, 0 + « 



Vobsjn, 0 
Vo 


(3) 


where n and / are the order and degree of the modes, respectively, 
b is the exponent calibrated from the Solar models, and a is a 
constant for the stellar model. The quantity Q„i presents the ratio 
between the inertia of the non-radial m ode /; and the iner tia of 
a radial mode Iq of the same frequency (lAerts et al.ll2010h . The 
inertia of mode 7/ is defined as; 


pR 

7, = 47r [Urf + 1(1 + l)Mrf]por^dr, (4) 

Jo 

where ^r(f) and ^h(r) are the radial and horizontal displacements 
of the modes at the radius r, respectively. The inertia ratio be¬ 
tween non-radial and radial modes can be expressed as 


Qni - IIIlo - 1 + 


\^h(r)\^Por^dr 

\&(rWpQr^dr 


(5) 


As discussed by lAerts et ^ (l2010h . the typical values for 
^h(R)l^r(R) are 0.01-0.1 for low-order p-modes, and 10-100 for 
high-order g-modes. This indicates that Qni ^ 1 for pure acoustic 
modes and Q„i > 1 forp-g mixed modes. 

To restrict the stellar parameters, we performed a Xc mini¬ 
mization by a comparison of models with observations; 


Tc 




C 


'theo 


■c° 


obs \ 


CT 


.obs 


(6) 


where C = (T’etr, [Fe/H], Av) and the cr*® denote the observa¬ 
tional errors. We chose models with < 1 as candidates for 
further pulsation analysis. In addition to the constraints from the 
atmospheric parameters and Av, we performed another xl min¬ 
imization by a comparison of the near-surface-corrected model 
frequencies with the observed frequencies; 


Av 


1 ^ 

--y 

N ^ 

l-l 


^theo _ yObs \ 


_obs 


(7) 


where the superscripts ‘obs’ and ‘theo’ correspond to individual 
frequencies from observations and modeling, respectively, and 
cr°'’* denotes the observational errors. We list the candidate mod¬ 
els v/ithxl^i, < 3000 for KIC 6442183 in Table|5l Since the / = 3 
modes of KIC 11137075 were not available, we chose the can¬ 
didate models with xl ^ 2000 for this star, as listed in Table 
|5] Note that the Brunt-Vaisala frequency affects the distributions 
of oscillation modes undergoing avoided crossings. Models with 
the similar global properties can have different mixed mode fre¬ 
quencies due to the diverse positions of avoided crossings. This 
results in the loxg^xlaii some models listed in Table |5] 

We note that the mass of these stars is not high enough to 
produce a convective core in the stellar interior, and the core 
overshooting will not have any impact on their structure and evo¬ 
lution. Therefore, we do not discuss the effects of overshooting 
in the following part for these low mass subgiants. 
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Table 5. The best candidate models for the two stars. The input parameters were; initial mass, mixing-length parameter (a), inclusion 
of diffusion (Diff) and initial metal abundance (Z,). The model parameters were: age (f), surface abundance ratio ((Z/Z)j), effective 
temperature (Teff), luminosity (L), radius (/?), gravity (logg), mean large separation(Av),;^^p,;^^^^^g and^y^n (columns 6-15). 


Model 

Mass 

(Me) 

a 

Diff 

z, 

t 

(Gyr) 

(ZIX), 

Teff 

(K) 

L 

(To) 

R 

{Ro) 

logg 

(dex) 

Ay 

(pHz) 

Xc 

Xtj=a 

^v,aU 

KIC 6442183 

1 

0.98 

1.70 

No 

0.012 

10.06 

0.017 

5703 

2.51 

1.63 

4.01 

64.86 

0.44 

29.85 

2514.27 

2 

0.99 

1.70 

No 

0.013 

9.95 

0.018 

5685 

2.50 

1.63 

4.01 

64.85 

0.27 

28.32 

599.80 

3 

1.00 

1.70 

No 

0.014 

9.80 

0.019 

5672 

2.49 

1.64 

4.01 

64.85 

0.34 

27.24 

77.20 

4 

1.04 

1.90 

No 

0.015 

8.62 

0.021 

5807 

2.82 

1.66 

4.01 

64.85 

0.49 

26.46 

70.98 

5 

1.03 

1.90 

Yes 

0.016 

8.85 

0.019 

5709 

2.62 

1.66 

4.01 

64.84 

0.09 

22.40 

1429.14 

6 

1.04 

1.90 

Yes 

0.015 

8.30 

0.017 

5783 

2.77 

1.66 

4.01 

64.84 

0.31 

27.61 

2208.02 

7 

1.04 

1.90 

Yes 

0.016 

8.50 

0.019 

5743 

2.69 

1.66 

4.01 

64.84 

0.03 

23.46 

532.24 

8 

1.04 

1.90 

Yes 

0.017 

8.68 

0.020 

5707 

2.63 

1.66 

4.01 

64.84 

0.15 

22.04 

46.23 

9 

1.04 

1.90 

Yes 

0.018 

8.85 

0.021 

5673 

2.57 

1.66 

4.01 

64.84 

0.57 

20.54 

571.86 

10 

1.05 

1.90 

Yes 

0.018 

8.50 

0.021 

5707 

2.65 

1.67 

4.02 

64.85 

0.32 

23.59 

1480.05 

KIC 11137075 

1 

0.96 

1.70 

No 

0.014 

11.60 

0.019 

5549 

2.19 

1.60 

4.01 

65.50 

0.26 

81.36 

1292.34 

2 

1.00 

1.90 

No 

0.015 

10.14 

0.021 

5686 

2.49 

1.63 

4.02 

65.51 

0.64 

64.19 

470.62 

3 

1.00 

1.90 

No 

0.016 

10.39 

0.022 

5645 

2.42 

1.63 

4.01 

65.52 

0.26 

43.86 

39.74 

4 

1.01 

1.90 

No 

0.017 

10.17 

0.024 

5645 

2.43 

1.63 

4.02 

65.53 

0.44 

19.19 

839.92 

5 

0.99 

1.90 

Yes 

0.016 

10.37 

0.019 

5587 

2.30 

1.62 

4.01 

65.53 

0.25 

19.16 

111.69 

6 

1.00 

1.90 

Yes 

0.016 

9.96 

0.019 

5622 

2.38 

1.63 

4.01 

65.48 

0.32 

135.53 

1712.46 

7 

1.00 

1.90 

Yes 

0.017 

10.16 

0.020 

5586 

2.32 

1.63 

4.01 

65.50 

0.05 

65.69 

431.52 

8 

1.00 

1.90 

Yes 

0.018 

10.36 

0.022 

5553 

2.26 

1.63 

4.01 

65.51 

0.10 

44.50 

32.11 


4.2. Optimal models 

The frequency separation between consecutive modes varies 
rapidly during an avoided crossing and so the mixed modes 
together with pure p-modes provide an important constraint to 
the stellar evolutionary s tate (IChristensen-Dalsgaard et aDl995t 
IPamvatnvkh et al.ll2004ll . We compared the frequencies and in¬ 
ertia ratios from observations and theoretical models listed in 
Table |5] The theoretical inertias for models were calculated us¬ 
ing equation!!] a nd the observed inertia ratio through the follow¬ 
ing function (e.g. lBenomar et aDl2014ll : 


curvature of the p-mode oscillation pattern. Although the dipole 
modes are p-g mixed modes in subgiants, we could estimate the 
expected pure acoustic modes v„^^i - v„p_o + 0.5Av-2Do through 
equations l9lor[T0l 

We define a quantity dvm^p to measure the frequency dif¬ 
ference between the mixed mode frequency v„_i and the pure 
acoustic modes Vn^.x. The quantity c/vm^p for dipole modes could 
be calculated using the following formula; 

= \Vn,X - (v«p,o + 0.5Av - 2Do)|. (11) 


h - y ^ /Te 

/o ‘AiVTi’ 


( 8 ) 


where A and T denote the amplitude and linewidth of modes, 
respectively, while the visibility Vi is the square root of the ratio 
between theoretical heights for the dipole and radial modes. We 
adopted the obs erved amplitude , linew idth and the value of Vi 
(y^=1.52) from iBenomar et ^ (1201 3h to deduce the observed 
inertia ratios. 

In solar-like stars, p- mode oscillatio ns are expected to follow 
the approximate relation (lTassoullll980l) : 

v„,; ~ (n + ^ + e)Av-1(1 H-l)Do, (9) 

where Dq is related to the interior structure of the star and the 
offset e is sensitive to the surface layers. The frequencies of the 
p-modes with the same degree show vertical ridges in the echelle 
diagram. The mean distance between dipole modes and radial 
modes with the same order would be 0.5Av - 2Dq, which is de¬ 
duced from th e approximate relation formulated in equation 0 
For red giants, iMosser et al.l (1201 Ibll introduced the expression 
for the pure p-mode eigenfrequency pattern: 


= («p + ^ + e)Av - 1(1 + l)Z)o H- ^(wp - nmax)^Av, (10) 


where is the p-mode radial order, Wmax is the order of the mode 
with highest power and a is a constant representing the mean 


The dvm-p is the frequency difference between the mixed modes 
and the nearest p-mode. For g-dominated mixed modes, the big¬ 
ger is the Qni, the stronger is the coupling between p- and g- 
modes. We used inertia ratios and <ivn,_p as criteria to constrain 
the stellar models of KIC 6442183 and KIC 11137075. 

Through the comparisons of theoretical and observed iner¬ 
tia ratios and dv^-p for KIC 6442183, we found that Models 3, 
4 without diffusion and Model 8 with diffusion match the ob¬ 
servations better than the other models. We estimated the fol¬ 
lowing parameters for KIC 6442183: M - 1.04;|;qq4Mo, R - 
1 ■ 66 ;|;qq 2 ^o and f = 8.65;(Qgg Gyr. By comparing Models 4 and 6 , 
we found that the Teff and luminosity estimated from the model 
without diffusion are higher than that estimated from diffusion 
models, which suggests that helium and heavy metal diffusions 
are not negligible. This is because the diffusion changes the 
opacity and equation of state for the models. The mass we find is 
approx imately 11% higher than that deduced bv iBenomar et al.l 
(l2Q14ll with the scaling relation. 

In Fig. 0 we show the results for KIC 11137075. As shown 
in the figure, the maximum of theoretical dipole t/vm-p is lower 
than the observations. For the I - 2 models presented in the 
third column of Fig. [ 8 ] Models 3, 5 without diffusion and Model 
8 with diffusion match the observations much better than the 
other four models, which means that the three models repro¬ 
duce the I - 2 avoided crossing better than the other four mod¬ 
els. Thus, Models 3, 5 and 8 were selected as the best models 
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Fig. 6. Comparison of echelle diagram, inertia ratios and 
dvm-p from observations and models without diffusion for KIC 
6442183. The color filled symbols denote the observed frequen¬ 
cies, and the unfilled symbols are the corrected theoretical fre¬ 
quencies. The squares denote the I - 2 modes, while the circles 
show / = 0, the triangle show I — 1 and the five-pointed star I - 3 
modes. We did not plot the error bars of frequencies since these 
are smaller than the symbol signs. 

for KIC 11137075. Finally, we obtained the following parame¬ 
ters for KIC 11137075; M = 1.00+° °|Mo, R = 1.63+° °|f?o and 
t - 10.36^q 2 o models with the mixing length a — 1.90 

reproduce the observations better than the models with a - 1.70. 

5. Discussions and Conciusions 

We carried out data processing and performed seismic analysis 
for subgiants KIC 6442183 and KIC 11137075, which were ob¬ 
served by the Kepler Mission. The main results and discussions 
are summarized as follows. 

We applied the Lomb-Scargle Periodogram to the corrected 
short cadence time series observed by Kepler to obtain the power 
spectra, and estimated the large frequency separation and fre¬ 
quency of maximum power for the two stars; 64.9 + 0.2 pHz and 
1225 + 17 pHz for KIC 6442183,65.5 ±0.2 ^uHz and 1171 +8 pHz 
for KIC 11137075, respectively. Individual mode frequencies 
were also extracted in this work. After carrying out asteroseismic 
analysis, we estimated the stellar parameters; M - 1.04^Q g 4 Mo, 
R = 1.66:““i/?o and t = 8 . 65!^-^2 Qyj. ^j. j^ic 6442183, 
M = 1.00:“;“jMo,= 1.63:“;°|/?o and f = 10.361“ Gyr for 



Fig. 7. Continuation of Fig. |6] 


KIC 11137075. Both the subgiants are shown to be solar-mass 
stars. 

For stars whose mass is similar to the Sun’s, the helium and 
heavy elements diffusion will affect the opacity and equation of 
state, which will change the thermal structure and affect the evo¬ 
lution of stars. Therefore, we have to consider the diffusion pro¬ 
cess in stellar models and seismic analysis. 

Pure p modes frequencies decre ase and g modes frequencies 
increa se with age for subgiants (IChristensen-Dalsgaard et all 
119951) . which results in the change of the frequency at which 
the avoided crossing occurs (lBedding|l2Q14l) . The frequencies of 
mixed modes sensitive to stellar interiors are useful to constrain 
stellar structures, especially for stellar interiors. The maximum 
of t/vm-p is coincident with the large inertia ratio, allowing us to 
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Fig. 8. Comparison of echelle diagram and dvm-p for observa¬ 
tions and modeling for KIC 11137075. The color filled symbols 
denote the observed frequencies, and the unfilled symbols are the 
corrected theoretical frequencies. The squares denote the / = 2 
modes, while the circles show I — 0, the triangles show I - \ 
modes. 


locate the mixed modes with the most g-mode property through 
these quantities. We hrst constrain stellar models with inertia ra¬ 
tios and dVm-p, which provides an impactful way to determine 
the parameters of evolved stars accurately. 

In future work, we plan to select some evolved stars with 
mixed modes such as subgiants, red giants and red clumps which 
constitute an evolution sequence. Then we compare the differ¬ 
ences of inertia ratios and t/vm-p between these stars and use 
these quantities to constrain the stellar interiors with the quanti- 
tis. 

Acknowledgements. We are grateful to the entire Kepler team. This work was 
supported by grants 10933002, 11273007 and 11273012 from the National 
Natural Science Foundation of China, and the Fundamental Research Funds for 
the Central Universities 


References 

Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, 
1st edn. (Springer) 

Aizenman, M., Smeyers, R, & Weigert, A. 1977, A&A, 58, 41 
Appourchaux, T., Gizon, L., & Rabello-Soares, M.-C. 1998a, A&AS, 132, 107 
Appourchaux, T., Rabello-Soares, M.-C., & Gizon, L. 1998b, A&AS, 132, 121 
Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 
Appourchaux, T., Benomar, O., Gruberbauer, M., et al. 2012, A&A, 537, A134 
Bahcall, J. N., Pinsonneault, M. H., & Wasserburg, G. J. 1995, Reviews of 
Modern Physics, 67, 781 

Barban, C., Deheuvels, S., Baudin, F., et al. 2009, A&A, 506, 51 
Basu, S., Turck-Chieze, S., Berthomieu, G., et al. 2000, ApJ, 535, 1078 
Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608 
Bedding, T. R. 2014, Asteroseismology, 22nd Canary Islands Winter School of 
Astrophysics, Cambridge University Press, Cambridge, p. 60 
Bi, S.-L., Basu, S., & Li, L.-H. 2008, ApJ, 673, 1093 
Bohm-Vitense, E. 1958, Z. Astrophys. 46, 108 

Borucki, W. J., Koch, D. G., Lissauer, J., et al. 2007, in Astronomical Society 
of the Pacific Conference Series, Vol. 366, Transiting Extrapolai* Planets 
Workshop, ed. C. Afonso, D. Weldrake, & T. Henning, 309 
Brandao, I. M., Dogan, G., Christensen-Dalsgaard, J., et al. 2011, A&A, 527, 
A37 

Brown, T. M., Gilliland, R. L., Noyes, R. W., & Ramsey, L. W. 1991, ApJ, 368, 
599 

Bruntt, H., Basu, S., Smalley, B., et al. 2012, MNRAS, 423, 122 
Benomar, O., Bedding, T. R., Stello, D., et al. 2012, ApJ, 745, L33 
Benomar, O., Bedding, T. R., Mosser, B., et al. 2013, ApJ, 767, 158 
Benomar, O., Belkacem, K., Bedding, T. R., et al. 2014, ApJ, 781, L29 
Chaplin, W. J., Houdek, G., Karoff, C., Elsworth, Y, & New, R. 2009, A&A, 
500, L21 

Chaplin, W. J., Appourchaux, T., Elsworth, Y, et al. 2010, ApJ, 713, L169 
Christensen-Dalsgaai'd, J., Bedding, T. R., & Kjeldsen, H. 1995, ApJ, 443, L29 
Christensen-Dalsgaard, J. 2004, Sol. Phys., 220, 137 
Datta, A., Mazumdar, A., Gupta, U., & Hekker, S. 2015, MNRAS, 447, 1935 
Demarque, R, Guenther, D. B., Li, L. H., Mazumdar, A., & Straka, C. W. 2008, 
Ap&SS, 316, 31 

Dotter, A., Chaboyer, B., Jevremovic, D., et al. 2008, ApJS, 178, 89 
Dupret, M.-A., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57 
Dziembowski, W. A., Gough, D. 0., Houdek, G., & Sienkiewicz, R. 2001, 
MNRAS, 328, 601 

Eggenberger, R, Charbonnel, C., Talon, S., et al. 2004, A&A, 417, 235 
Eerguson, J. W., Alexander, D. R., Allard, R, et al. 2005, ApJ, 623, 585 
Garcia, R. A., Mathur, S., Salabert, D., et al. 2010, Science, 329, 1032 
Gai'ci'a, R. A., Hekker, S., Stello, D., et al. 2011, MNRAS, 414, L6 
Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010a, PASP, 
122, 131 

Gilliland, R. L., Jenkins, J. M., Borucki, W. J., et al. 2010b, ApJ, 713, L160 
Guenther, D. B. 1994, ApJ, 422, 400 

Huber, D., Stello, D., Bedding, T. R., et al. 2009, Communications in 
Asteroseismology, 160, 74 
Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 
Kallinger, T., Weiss, W. W., Barban, C., et al. 2010, A&A, 509, A77 
Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87 

Kjeldsen, H., Bedding, T. R., & Christensen-Dalsgaard, J. 2008, ApJ, 683, L175 
Kjeldsen, H., Christensen-Dalsgaard, J., Handberg, R., et al. 2010, 
Astronomische Nachrichten, 331, 966 
Lomb, N. R. 1976, Ap&SS, 39, 447 






















Tian et al.: Asteroseismic analysis of subgiants KIC 6442183 and KIC 11137075 


Markwardt, C. B. 2009, Astronomical Data Analysis Software and Systems 
XVIII, 411,251 

Markwardt, C. 2012, Astrophysics Source Code Library, 1208.019 
Mathur, S., Handberg, R., Campante, T. L., et al. 2011, ApJ, 733, 95 
Metcalfe, T. S., Monteiro, M. J. R F. G., Thompson, M. 1., et al. 2010, ApJ, 723, 
1583 

Molenda-Zakowicz, J., Sousa, S. G., Frasca, A., et al. 2013, MNRAS, 434, 1422 

Montalban, J., Miglio, A., Noels, A., et al. 2013, ApJ, 766, 118 

Mosser, B., & Appourchaux, T. 2009, A&A, 508, 877 

Mosser, B., Barban, C., Montalban, J., et al. 201 la, A&A, 532, A86 

Mosser, B., Belkacem, K., Goupil, M. J., et al. 201 lb, A&A, 525, L9 

Mosser, B., Benomar, O., Belkacem, K., et al. 2014, A&A, 572, L5 

Osaki, J. 1975, PASJ, 27, 237 

Regulo, C., & Roca Cortes, T. 2007, A&A, 469, 233 

Pamyatnykh, A. A., Handler, G., & Dziembowski, W. A. 2004, MNRAS, 350, 

1022 

Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 

Roxburgh, I. W. 2009, A&A, 506, 435 

Roxburgh, I. W., & Vorontsov, S. V. 2006, MNRAS, 369, 1491 

Scargle, J. D. 1982, ApJ, 263, 835 

Tassoul, M. 1980, ApJS, 43, 469 

Thompson, B., Frinchaboy, R, Kinemuchi, K., Sarajedini, A., & Cohen, R. 2014, 
AJ, 148, 85 

Thoul, A. A., Bahcall, J. N., & Loeb, A. 1994, ApJ, 421, 828 
Tian, Z. J., Bi, S. L., Yang, W. M., et al. 2014, MNRAS, 445, 2999 


9 



